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Abstract 

We combine the fixed-order evaluation of the bb sum rules with a non- 
relativistic effective-theory approach. The combined result for the n-th 
moment includes all terms suppressed with respect to the leading-order 
result by 0(a? s ) and O [^a s y/n) l a^) , counting a s ^/n ~ 1. When com- 
pared to experimental data, the moments thus obtained show a remark- 
able consistency and allow for an analysis in the whole range 1 < n < 16. 



1 Introduction and outline 



Near threshold, the cross section for the production of a bb pair, a(e + e~ —>■ bb), 
is extremely sensitive to the mass of the bottom quark m&, which allows for a 
precise determination of m&. This is usually done by considering sum rules [1] 
and defining the n-th moment 

dS T, , S ( d V TT/ 2M 



where II(g 2 ) is the vacuum polarization, e& = —1/3 the electric charge of the bot- 
tom quark and Rbh(s) = cr(e + e~ — >■ bb)/cr(e + e~ — > n + ^ ) the normalized cross 
section. In order to extract m&, the theoretical evaluation of M n is compared 
to the experimental value. From the experimental point of view, the moment 
obtains contributions from the six T bound states and from the continuum cross 
section above threshold. For increasing n, the contribution from the experimen- 
tally poorly known continuum cross section becomes less and less relevant due 
to the suppression l/s n+1 . As for the choice of the parameter n, there are two 
complementary approaches. Either n is assumed to be rather small, i.e. n < 4 in 
which case II (q 2 ) is computed in a standard weak coupling fixed-order approach 
or n is assumed to be rather large n > 8 in which case the moments are evaluated 
in a non-relativistic effective-theory approach. 

In the standard fixed-order (FO) approach, the vacuum polarization is written 

as 

where N c = 3 is the colour factor and the coefficients C n are evaluated as a series 
in the strong coupling a s . These coefficients depend on the mass scheme that is 
used. We will indicate this dependence by a label X, i.e. m x denotes the bottom 
quark mass in a particular scheme and C ny x are the corresponding coefficients. 
From the knowledge of C n ^x, the moments are obtained as 

M n x = ~N C ei - — -T5- C n x (3) 
' 4 b (2m x ) 2n 

Being observables, the moments should be scheme independent. However, since 
the perturbative series is truncated, there is a residual scheme dependence left 
in M n which again is indicated by the label X. The coefficients C n have been 
computed up to 0(a 2 ), i.e. three loops for n < 8 in Ref. [2] and up to n < 30 
in Ref. [3]. The four-loop coefficient is known for n — 0,1 [HE] and these results 
have been used to obtain precise values for m = the bottom quark mass in 
the MS-scheme [6] . 
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At I + 1 loops, the coefficients C n contain terms n~^' 2 (a s y/n) 1 . Thus, if n 
increases, the higher-order terms become more important and for ^fri ~ a~ l the 
standard fixed-order approach completely fails. This is related to the fact that 
in a strict expansion in a s the theory does not contain bound states. Given that 
for increasing n the moments are dominated by the lowest resonances it is thus 
not surprising that a FO approach does not very well describe M n for large n. 
As a consequence, mass determinations using this approach [7J [6] use small 
values of n. 

In order to describe the weak coupling bound states in the bb system, we 
have to consider the non-relativistic sum rule. The starting point is the solu- 
tion to the Schrodinger equation describing a non-relativistic bb pair interacting 
through the Coulomb potential — CVa s /r with Cf = 4/3 a colour factor. This 
resums all terms of the form v (a s /v) 1 in R b i, where v is the small velocity of the 
heavy quarks and, therefore, resums all terms of the form n _3//2 (a s y/n) 1 in M n . 
Higher order corrections in a s as well as v are taken into account using Quantum 
Mechanics perturbation theory and counting a s ~ v. This is done most effi- 
ciently in the framework of an effective theory (for a review see Ref. [9]). Within 
the effective-theory (ET) approach the next-to-next-to-leading order (NNLO) 
corrections to the non-relativistic sum rules, including all terms suppressed by 
a 1 ~ a s/ V™ ~ n ~ l with respect to the leading-order result, have been computed 
and used to determine the bottom quark mass [10]. The theoretical predic- 
tions can be improved upon by resummation of large logarithms [llj of the form 
{a s \ogv) [ and including these terms in the non-relativistic sum rule leads to a 
much more robust determination of the bottom quark mass [T2] . 

Ultimately the most interesting quantity is fn, the bottom quark mass in the 
MS-scheme. Using non-relativistic sum rules, fn has to be determined in two 
steps. First, a so called threshold mass p21[TlJ[T5] has to be used which is closely 
related to the pole mass but accounts for the cancellation of the renormalon 
ambiguity in the observable M n . In a second step, the threshold mass is related 
to fn. Using this approach, the moments can be determined reliably for large 
values of n, as long as non-perturbative contributions are not too important. 
However, for small values of n this approach breaks down due to the neglect of 
terms suppressed for large n or small v. As an example consider terms of order 
atgV 3 in R b i which are kept at NNLO in an ET approach, whereas terms of order 
a.gV 5 are dropped. These terms result in contributions of the order a°n~ 5 ^ 2 and 
a o n -7/2 reS p ec tively in M n . It is clear that the latter are suppressed in the non- 
relativistic sum rule, i.e. for large n, but their neglect invalidates the n — > 1 limit 
of the result obtained in the non-relativistic approach. 

Large n and small n applications of the sum rules and the corresponding 
determinations of m& both have their advantages and disadvantages. From the 
large n point of view, one advantage is that due to M n ~ l/{2mb) 2n the moments 
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are much more sensitive to the bottom quark mass for large n. Also, they are 
virtually insensitive to the continuum contribution. Since this contribution is 
experimentally only known very poorly, small n determinations of the bottom 
quark mass crucially rely on the precise treatment of the data in the threshold 
region and have to use perturbative QCD input for R h i above threshold. Given 
that the experimental error in the m& determination for small n is dominant, 
rather subtle changes in the treatment can have significant effects on the extracted 
value of m;, and, in particular, its error. On the other hand, the perturbative series 
is much better behaved for small n. In fact, the non-relativistic sum rules suffers 
from very large corrections and even though the resummation of the logarithms 
substantially improves the behaviour of the perturbative series, the situation is 
far from ideal, and the dominant error still comes from the neglect of higher-order 
corrections. For completeness we mention again that in the large n approach m 
cannot be obtained directly but only through the intermediate use of a threshold 
mass. However, this is not necessarily a disadvantage, since threshold masses are 
important and useful in their own right. 

It is natural to ask whether it is not preferable to combine both approaches 
and perform an all n analysis of the sum rule. Since the two approaches use 
different techniques and are sensitive to rather different experimental input it 
certainly would give increased credibility if the extracted value of m& varies very 
little with n, and it would allow to get a better handle on the determination of 
its error. The choice of n in such a combined analysis is only limited by the 
non-perturbative corrections. In order to get an estimate of the importance of 
these corrections it is useful to consider the contribution of the gluon condensate 
to the sum rule [161 EJ- Even though this contribution grows rapidly with n, for 
realistic values of the gluon condesate it is below 0.1% for n < 12 and reaches 
about 1% for n = 16. It is thus legitimate to neglect non-perturbative corrections 
to the 6-quark sum rules as long as n is not chosen to be too large. 

In this paper we consider the first 16 moments, starting in Section [2] with the 
fixed-order approach. Even though the large-n behaviour of the FO results is 
better than anticipated, will find the expected problems for large n and turn in 
Section [3] to the non-relativistic sum rule in order to illustrate its behaviour as a 
function of n. Finally, in Section HI we combine the two approaches by adding to 
the non-relativistic sum rules all terms order that have been missed. As we 
will see this allows to obtain a consistent description of M n for 1 < n < 16. We 
refrain from presenting another extraction of rat,, since all theory input used in 
this analysis has already been used for a bottom mass determination [51 [6j [12]. 
The main aim of the paper is to establish the fact that a future analysis, once 
further improved theoretical results and hopefully better experimental data is 
available, should consider the full range of n in order to get a better control of 
the different systematic uncertainties. 
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2 Fixed order results, small n 



In a fixed-order approach, M n x can be written as in Eq. 
C n x has the structure 



and the coefficient 
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In the on-shell scheme we use the notation Cn^ = C^os 



and we have 



(ii) 



<*J> with / > 1, 



Cn 22 ^ = Cn 33 '* = 0. The logarithmic coefficients of order a k s , i.e. C, 
can be predicted from the lower order coefficients of order a™, m < k. The 
three-loop coefficients C^x have been computed up to n = 8 in Ref. |2j and later 

up to n = 30 in Ref. [5]. The four-loop coefficient C[ 3 X is also known [U E] but 
for n > 1 these coefficients have not yet been computed. 

The relation between C^ x and Cn l \ the coefficients in the scheme X and 
the on-shell scheme respectively, where the corresponding masses (we denote the 
pole mass by m = mos) are related by 
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It is clear that theses relations break down for large n, due to the terms of 



order (na* 



In fact, the behaviour for n — > oo seems to be even worse than 
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n~ z l 2 (a s yjn) 1 as mentioned in the introduction. This is due to the shift in the 
mass scheme. To be precise, C n , the coefficients in the on-shell scheme behave 
like n~^l 2 (a s ^Jn) l for n — > oo. For any threshold scheme X, the factors 5m^ 
have to have an additional suppression Sm^ ~ oi s in order not to destroy the 
behaviour n~ 3 / 2 (a s y/n) 1 for n — > oo. For the MS-scheme this is not the case and, 
as we will see, this scheme is particularly inappropriate for large n. 
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Figure 1: Scale dependence of (some of) the first 16 moments in the MS- 
scheme. The moments are evaluated with m = 4.184 GeV in a fixed-order 
approach including terms up to 0{a 3 ). 




In order to substantiate this point consider the scale dependence of the first 
16 moments evaluated in a fixed-order approach up to 0(a 3 ) in the MS-scheme. 
The n-th moment has mass dimensions [m^] -2 ™ and the moments in this paper 
are always given in units [GeV]~ 2n . Since C^ 3 ^ is not known for n > 1 we set 

C^ 3 ^ = C| 3 ^- We fix m = ^Ms( m Ms) = 4.184 GeV and then evaluate ^ms(a*) 
using the renormalization-group equations to four- loop accuracy [18]. We then 
use this value for the mass and the scale \i to evaluate M n ^ and vary the scale 
in the region 4 GeV < \i < 10 GeV. As can be seen in Figure [T], for n < 6 the 
moments are very stable for the whole range of /i, but for larger values of n the 
scale dependence deteriorates rapidly and it is clear that for n > 8 no information 
can be extracted from these results any longer. The precise shape of the curves in 
Figure [T] depends to some extent on details as how to treat the flavour threshold 
and the value of C*- 3 ^, but the main point is not affected by these issues. 

The same exercise can be repeated for a threshold mass. As discussed above in 
this case we would expect a somewhat better behaviour for large n. To investigate 
this, we use the PS-scheme [H] and set mps = mpg(/i^ = 2 GeV) = 4.505 GeV 
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which, using three-loop conversion, corresponds to m — 4.184 GeV. Since Srripg = 
Cf IXf/^ps we have to choose the factorization scale ilf ~ rriba s to ensure the 
additional suppression 5m^ ~ a s mentioned above and the standard choice is 
ix p = 2 GeV. Again we evaluate the first 16 moments varying the scale in 
the region 2.5 GeV < ll < 20 GeV. The results depicted in Figure [2] show a 
remarkable stability with respect to the scale variation. The scale dependence 
does increase for increasing n but remains in much better control than in the 
MS-scheme. 




H [GeV] 

Figure 2: Scale dependence of (some of) the first 16 moments in the PS- 
scheme with mps = 4.505 GeV. The moments are evaluated in a fixed- 
order approach including terms up to 0{o? s ). 

Even though the results shown in Figures [1] and [2] nicely confirm our expec- 
tations it is clear that even in the PS-scheme the results become unreliable for 
large n. To obtain a more complete picture we now turn to the evaluation of 
the moments in the non-relativistic effective-theory approach and compare these 
results with the fixed-order results of this section. 

3 Effective theory results, large n 

As mentioned in the introduction, in the effective-theory approach we start from 
the Schrodinger equation describing a non-relativistic heavy quark pair with en- 
ergy E = \fs — 2ml) interacting through the potential — CVa s /r. The cross section 
R{s) is related to the imaginary part of the corresponding Green function at the 
origin. Working in dimensional regularization in d = 4 — 2e dimensions and 
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minimally subtracting the 1/e ultraviolet singularity, the leading order Coulomb 
Green function at the origin is given by [19] 



Gf (0, 0; E) = (- + - log _ _ + 1E + , (1 _ A) j 

where A = Cp ot s / (2 a/— Ejm^) and the leading-order cross section is given by 

e 2 

#(£) = 67riV c -\ Im [G c (0, 0; E)\ (12) 

Higher-order corrections are computed by perturbative insertions of higher-order 
corrections to the potential. The moments are then evaluated preforming the 
integration indicated in Eq. (pQ). In the literature different options on how to 
treat the prefactor l/s n+l have been used. Either, this factor can be expanded, 
writing 

f°° 2dE / E nE 2 A . . 

M n = / — - - e m b 1 + + . . . R(E) 13 

J_ 00 (2m b y^ V 2m b (2m 6 ) 2 ) K ' K ' 

where the ellipses stand for higher-order terms in the non-relativistic expansion, 
or it can be left unexpanded 

/oo r> j 771 

Eqs. f|T3|) and f|T4l) agree at NNLO in the effective theory, but will differ consid- 
erably for small n. In this section we will use the strictly expanded approach, 
Eq. f|T3|) . but we will come back to this issue in Section HI 

In Figure [3] we show the scale dependence of (some of) the first 16 mo- 
ments evaluated in the effective theory with mps = 4.505 GeV. These results 
are complete at NNLO. Thus, counting a s ~ they include all terms scal- 

ing like n~ 3 / 2 (a s ^/n) 1 a 2 . Furthermore large logarithms are resummed counting 
a s logn ~ 1. The results presented here are complete at next-to- leading logarith- 
mic (NLL) accuracy and contain some known contributions at next-to-next-to- 
leading logarithmic (NNLL) accuracy [12j. The scale dependence of the first few 
moments is very strong, indicating the expected breakdown of the ET approach 
for small n. However, for n > 5 the results are very stable. Increasing n further 
to n > 14 leads to an enhanced scale dependence. This is not unexpected, since 
missing higher-order and NNLL terms as well as non-perturbative corrections 
become increasingly important. 

Of course, the scale dependence is at best a very rough indicator of the re- 
liability of the results. In the following section we will combine the results of 
Sections [2] and [31 including all terms 0(a^) of the FO approach and all (known) 
NNLL terms of the ET approach and investigate the relative importance of the 
various corrections. 
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Figure 3: Scale dependence of (some of) the first 16 moments in the 
PS-scheme with mps = 4.505 GeV. The moments are evaluated in an 
effective-theory approach including terms up to NNLL accuracy. 



4 Combined analysis 

In this section we present the results of a combined approach, i.e. results that 
are complete at 0(a^) in the FO approach and complete at NNLL in the ET 
approach. The results are obtained simply by adding the FO and ET results and 
subtracting the doubly counted terms. In order to obtain the doubly counted 
terms we expand the ET result in the coupling a s and retain all terms of 0(a^). 
Note that these terms depend on the precise implementation of the non-relativistic 
expansion. In particular, they depend on whether Eq. f|T3|) or Eq. (fT4l) is used. 
Any implementation that is equivalent at large n of the ET result can be used, 
as long as the subtraction terms are treated consistently. 

Using the implementation according to Eq. ( {TBI and expanding the ET result 
in a s leads to integrals of the form 

f°° 2dE _^ \fm b \ x , k f-Em b \] 

for the corresponding contribution to M n , where E = E + i0 + is understood. 
These integrals can either be computed numerically or obtained analytically by 
differentiation with respect to y of 
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The corresponding integral using the implementation according to Eq. (j!4p is 
given by 
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In the derivation of these results we assumed that \i is independent of E. 



(17) 



n 


1 


2 


3 


4 


5 


6 


7 


'exact LO' 


9.235 


5.021 


3.998 


3.701 


3.713 


3.914 


4.267 


0(a° s ) 


5.458 


2.377 


1.594 


1.275 


1.124 


1.053 


1.030 


0{a\) 


8.294 


4.124 


3.028 


2.600 


2.430 


2.394 


2.445 


0{al) 


9.065 


4.795 


3.704 


3.321 


3.224 


3.287 


3.463 


0{al) 


9.211 


4.976 


3.926 


3.595 


3.561 


3.702 


3.975 


0{<4) 


9.233 


5.014 


3.983 


3.676 


3.673 


3.853 


4.176 


0(a 5 s ) 


9.236 


5.020 


3.995 


3.696 


3.704 


3.899 


4.242 



Table 1: Comparison of the exact LO result in the effective theory to the 
expanded results, Eq. (Tl9j) . All entries are multiplied by 10 2n+1 . 



In order to illustrate the procedure, let us take the leading-order result in the 
effective theory, given in Eqs. (TTIi) and ( TT2l) and expand it in a s to say 0(a b s ) 



R(E) = ^N c C F e 



3 of 1 ( , -4m b E\ 



j | -]o R — -^1 + 2a 2 J^ 1 \l) (18) 



«^V 2) (i) + l -a±e v (3) (i) - ^*V 4) (i) + • • • 



where we introduced I = X/a s = CV/(2a/— E/m b ) and ip^ denotes the k-ih 
derivative of the ^-function. Integrating this series according to Eq. (fl3|) . drop- 
ping the higher-order terms in E/m b , we get 



3n" 3 / 2 



4(2m 6 ) 



2n ' 



(19) 



- « 3 \ ^ (2) (1) + « 4 ^ ^ (3) (1) - « 5 ^ ^ (4) (1) + • • ■ 
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with a = Cir(a sV /n). We can now check how the expanded result, Eq. (fi~9l) 
approaches the 'exact' leading-order result in the effective theory. This is done 
in Table [U where we show the results of performing according to Eq. (fl3l) the 
integration of R{E) as given in Eq. (Tl2l) . As in the derivation of Eq. (TT9l we drop 
higher-order terms in E/rrib and, for convenience, multiply by 10 2n+1 . The results 
for n < 7 with rrib = 4.505 GeV and [i = 4.5 GeV are shown in the second row, 
labelled 'exact LO'. The other rows contain the successive approximations given 
in Eq. ffT9l . with a s = a s (/i = 4.5 GeV) = 0.2198. As expected, the expanded 
results approach the 'exact result' faster for smaller values on n. Including all 
terms up to 0(a 3 s ) the relative error is {0.3%, 0.9%, 1.8%, 3.0%, 4.3%, 5.7%, 7.4%} 
for n = {1, 2, 3, 4, 5, 6, 7} respectively. The doubly counted terms to be subtracted 
in the combined analysis at this order correspond to the terms given in the row 
labelled 0(a 3 8 ). 
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Figure 4: Scale dependence of the first 16 moments in the PS-scheme with 
mps = 4.505 GeV. The moments are evaluated in a combined approach 
including terms up to NNLL accuracy and order O(a^). 

Repeating this exercise with the full NNLL effective-theory result and com- 
bining the FO and ET results using the definition Eq. ([TBI we evaluate again 
the first 16 moments and depict their scale dependence in Figure HI Comparing 
Figures H] and [3] we note that, as expected, the difference is small for large n and 
large for small n. The FO corrections to the ET results are < 5% for n — 10 
(except for very small scales), increasing to < 10% for n = 4. For n < 2 the 
corrections completely change the shape of the curve, indicating the importance 
of relativistic corrections to the non-relativistic sum rules. On the other hand, 
comparing Figures H] and [2] the situation is just reversed. For small n the correc- 
tions are small (except for very small scales) and they increase with increasing 
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n 


1 


2 


3 


4 


6 


8 


10 


13 


16 


10' 2n+1 M n 


4.70 


2.94 


2.42 


2.23 


2.21 


2.42 


2.79 


3.68 


5.13 


comb/ET 


1.61 


1.03 


0.98 


0.97 


0.97 


0.97 


0.97 


0.97 


0.99 


comb/FO 


1.05 


1.07 


1.09 


1.10 


1.14 


1.16 


1.19 


1.22 


1.24 



Table 2: Comparison of the combined evaluation of selected moments 
M n with the ET and FO approach. The moments are evaluated with 
mps = 4.505 GeV with the scale (i = 4.5 GeV. 



n indicating the importance of resumming terms (a 3 y/n) for large n. This is 
also confirmed by Table El where we list the value of some combined moments 
(second row), as well as the ratio of the combined moment to the ET result (third 
row) and FO result (fourth row) respectively. We should stress that the ratios 
in Table [2] depend on the scale choice \i — 4.5 GeV and only give an incom- 
plete picture. In particular, the corrections to the ET result for n G {2,3,4} are 
larger than what might be inferred from Table El This is illustrated in Figure [5] 
where the scale dependence of the second moment is plotted and compared to 
the experimental moment with its error, indicated by the black line and the grey 
rectangle. The FO result is plotted as the light blue line. The ET result is eval- 
uated using Eq. ( fi~3l (solid dark blue line) and Eq. ( fT4l) (dashed dark blue line). 
As mentioned above, the two implementations differ considerably (for small n). 
However, the corresponding combined results, depicted as solid and dashed ma- 
genta lines respectively, are virtually independent of the implementation, since 
differences in treating higher-order in n terms are compensated for by adding the 
full n dependence up to 0{pP s ) through the FO result. From Figure [5] we can also 
see that the value of 1.03 given in Table El for the ratio of the combined and ET 
result for the second moment is a coincidence of the scale choice \l — 4.5 GeV 
and not necessarily indicative of the typical size of the corrections. 

The experimental moments used in Figure and the following plots have 
been determined by taking into account M™ s , the contribution due to the six 
lowest resonances, and using perturbative QCD in the region > 11.2 GeV to 
obtain the continuum contribution M™ nt . This follows closely Ref. [7] from which 
we also adopt the treatment of M^ n , the additional contribution in the region 
y's > 11.2 GeV. Due to the uncertainty and the lack of precise experimental 
data in the region around and just above threshold we add the errors linearly. 
The experimental moments and their errors are listed in Table |3j The first four 
moments agree within errors with those given in Ref. [6], but our experimental 
moments have a larger error. 

In order to get a better understanding of the relative importance of the var- 
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effective theory 
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fixed order 
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fi [GeV] 

Figure 5: Scale dependence of the second moment in the PS-scheme eval- 
uated using a FO (light blue curve) , ET (magenta curves) and a combined 
approach (dark blue curves). The dashed curves have been obtained using 
Eq. (fTB|) . whereas the solid curves have been obtained using Eq. (fl4j) . 



n 


1 


2 


3 


4 


6 


8 


10 


13 


16 


10 2n+l M e x p 


4.51 


2.81 


2.31 


2.13 


2.11 


2.30 


2.64 


3.40 


4.53 




0.15 


0.12 


0.10 


0.08 


0.07 


0.06 


0.06 


0.06 


0.07 



Table 3: Values of some selected experimental moments and their errors. 



ious corrections and the range of applicability of the various approximations we 
compare the first 16 moments in the PS-scheme to the experimental moments. 
We evaluate the moments in the FO, ET and combined approach as well as in 
the MS fixed-order approach. Note that m = 4.184 GeV has been determined 
by requiring the first moment in the MS fixed-order approach to agree with the 
experimental value. This then fixes mps = 4.505 GeV and all further moments. 
For the PS-scheme we vary the scale in the range 2.5 GeV < \i < 10 GeV, 
whereas for the MS scheme we vary the scale in the range 4 GeV < // < 10 GeV 
as explained in Section [TJ The results are shown in Figure El Note that the scale 
fi — 4.5 GeV chosen in Table [2] is close to the upper end of the scale variation 
shown in Figure [61 This explains why the values in Table [2] are generically larger 
than the experimental moments given in Table [3j The most striking feature of 
Figure [6] is that in the PS-scheme all three approaches give very similar results. 
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In particular, the FO approach gives good results even for n = 16. The ET 
approach seems to be valid down to n = 3 and only breaks down for n < 2. We 
note this seems to be a general feature of any suitably defined threshold mass. In 
particular, we have checked that for the RS-mass [15] the results are very similar. 
On the other hand, the situation is rather different in the MS scheme. The FO 
approach gives excellent results for n < 7 and then breaks down abruptly. This 
can also be inferred from Figure [H Had we chosen to limit the scale variation by 
say \x < 7 GeV we would have obtained good results up to n = 9. 



PS fixed order □ experiment 
PS effective theory 
PS combined 
MS fixed order 
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Figure 6: Comparison of the experimental moments to the PS-scheme 
calculation in a FO approach (left /light blue bands), in an ET approach 
(middle/magenta bands) and in a combined approach (right/dark blue 
bands). The bands have been obtained by varying 2.5 GeV < fj, < 10 GeV. 
Similar bands for the MS-scheme are shown in red. 



Finally, we present a similar plot for the OS scheme. We fix the value of 
the pole mass to make the first moment in the FO approach to agree with the 
experimental moment. This results in m = 4.85 GeV. The we proceed as in the 
case of the PS scheme. As mentioned at the beginning, the OS scheme is not 
well suited for a precise determination of quark masses and we would expect the 
results to be less consistent than with other mass definitions. This is what we 
find in Figure [7J The FO results are inconsistent with the experimental values of 
the moments for n > 6. Accordingly, the corrections (a s y/n) 1 are more important 
than in the PS scheme and bring the combined results into agreement with the ET 
results and the experimental values. The ET results agree with the experimental 
moments for all values of n, but for n = 1 the scale dependence is enormous, 
making the result meaningless. Overall, the scale dependence is considerably 
larger than in the PS scheme, in agreement with our expectations. 
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Figure 7: Comparison of the experimental moments to the moments com- 
puted in the OS scheme with the pole mass set to m = 4.85 GeV in a FO 
approach (left/light blue bands), in an ET approach (middle/magenta 
bands) and in a combined approach (right/dark blue bands). The bands 
have been obtained by varying 2.5 GeV < fj, < 10 GeV. Similar bands for 
the MS-scheme are shown in red. 



5 Conclusions and outlook 



The main result of this analysis is that there is no need to make the standard 
separation into large-n and small-n analyses of the bb sum rules. For a suitably 
defined threshold mass, the fixed-order results are remarkably consistent even for 
large values of n. With hindsight one might argue that from a numerical point 
of view the expected breakdown does not happen at a s ^fn ~ 1 but rather at 
(a s /ir) ^fn ~ 1. Even in the case of the MS-mass where the large-n behaviour is 
worse due to the presence of terms (a s n)\ a fixed-order approach is applicable for 
values of n up to n ~ 6. On the other hand the non-relativistic sum rule can also 
be applied for values of n that are much smaller than what naively could have been 
expected. Overall, we obtain a very consistent picture. With the availability of 
the NNNLO corrections in the effective theory [20] and the prospect of complete 
results in the fixed-order approach at 0(aP s ) also for n > 1, the sum rule is likely 
to be the observable of choice for bottom quark mass determinations. Due to the 
large range of n that can be used, non-perturbative corrections are well under 
control and additional effects such as non-vanishing charm mass [2TJ [8] can be 
included as well. In view of the progress on the theoretical side, more precise 
experimental data of the 6-quark cross section just above threshold would be 
most welcome. 
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